Image-potential induced states at metal surfaces 
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The unequivocal observation in the early 80's of a new class of unoccupied surface states that are 
bound to the vacuum level of a variety of metal surfaces has led over the last two decades to an 
active area of research in condensed-matter and surface physics. Here we discuss the key ingredients 
of the theory of these so-called image states, which are originated in the long range nature of the 
image potential outside a solid. 
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I. INTRODUCTION 

It was shown by Cole and Cohen Q and by Shikin 
H that for negative electron affinity materials the im- 
age potential outside the surface should trap electrons in 
Rydbergs states, thereby explaining the experimentally 
observed trapping of electrons at the surface of liquid he- 
lium H 4 , J| . These so-called image states were then 
invoked [3, uj to explain the fine structure that had been 
observed in the difraction of low-energy electrons from 
metal surfaces with a band gap near the vacuum level 
H p|, p^|. Echenique and Pendry [jll) also investigated 
the observability of image states via low-energy-electron 
difraction (LEED) experiments, and a discussion of the 
lifetime broadening of these states led them to the impor- 
tant conclusion that they could in principle be resolved 
for all members of the Rydberg series. 

Another type of surface states in the band gaps of free- 
electron-like s,p bands, which can occur even for a step 
barrier in the absence of the image potential, were pre- 
dicted by Shockley j^] . These states have been observed 
with the use of photoemission in a variety of metal sur- 
faces 0,0. 

In angle-resolved photoemission experiments a sample 
is illuminated with monochromatised photons and the en- 
ergy spectrum of the photoemitted electrons is examined, 
thus providing information about the density of occupied 
states. In order to investigate unoccupied states in the 
range between the Fermi and the vacuum levels, one has 
to use other methods such as inverse photoemission [[18] . 
In these experiments electrons with a given energy are 
incident on a solid surface and the energy of the emit- 
ted photon is measured. Hence, this technique allows to 
map out the features of the unoccupied density of states 
and to derive the energy and momentum of bound states 
by simply measuring the energy and momentum of the 
incident electron and the energy of the emitted photon. 

Johnson and Smith ]l9| pointed out that image- 
potential induced bound states were potentially observ- 
able by angle-resolved inverse photoemission. By us- 
ing this technique, Dose et al. [p0| and Straub and 



Himpsel |2l| reported the first conclusive experimental 
evidence for image-potential bound states at the (100) 
surfaces of copper and gold. Since then, several ob- 
servations of image states were made using this tech- 
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|24| , |25| , |26| , |27| , |28| , and also the new 
high-resolution techniques of two-photon photoemission 
(TPPE) HI |(], H and time-resolved two-photon pho- 
toemission (TR-TPPE) @ |§. In TPPE, intense laser 
radiation is used to populate an unoccupied state with 
the first photon and to photoionize from this interme- 
diate state with the second photon. In TR-TPPE, the 
probe pulse which ionizes the intermediate state is de- 
layed with respect to the pump pulse which populates 
it, thus providing a direct measurement of the intermedi- 
ate state lifetime @ [|| |^ §7|, §§, |H @- Since image 
states are typically decoupled from bulk states, their long 
lifetimes make them to be ideally suited as intermediate 
states. Both image states and the image-like form of the 
surface barrier outside a solid also play a role in scanning 
tunneling microscopy (STM) @ || || ||. 

The binding energies of image sates were originally 
investigated in terms of a simple multiple-scattering 
approach in which surface states are viewed as waves 
that are trapped and multiply scattered between the bulk 
crystal on one side and the surface barrier on the other 
Jio| , pH . Soon after the first observation of these states, 
it was shown that the crystal-face dependent binding en- 
ergies of both Shockley and image states could be under- 
stood in terms of the multiple-scattering approach [[f6| , 
and this approach was then combined with a nearly-free- 
electron (NFE) description of the bulk band structure 
i pTf to predict Shockley and image-state binding energies 
that were in reasonable agreement with photoemission 
and inverse photoemission measurements KM |4^, ^l) . 
Wave-function matching techniques were also applied to 
obtain both binding energies and wave functions of sur- 
face states and to demonstrate that the spatial extent of 
these states is determined by their binding energy |5^] . 

Since the pioneering work of Appelbaum and Hamann 
p3j , first-principles calculations of the electronic struc- 
ture of solid surfaces have been carried out within the 
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framework of density-functional theory (DFT) pi] by 
replacing the full many-body Schrodinger equation by 
an effective one-particle equation that is typically solved 
sclf-consistently in the so-called local-density approxima- 
tion (LDA) |35|]. Although this approximation ceases 
to be valid in the tail of the electron density where 
the asymptotic value of the one-particle LDA effective 
potential does not reproduce the correct image-like be- 
haviour, LDA calculations have proved to be succesful 
in the description of both bulk states travelling away 
from the surface and Shockley surface states decaying 
exponentially away from the surface. Nevertheless, im- 
age states can only be described within a correct non- 
local description of the one-particle potential in the vac- 
uum side of the surface. By matching a long-range im- 
age potential onto the self-consistent LDA crystal po- 
tential, the experimentally observed binding energies 
and effective masses of image states on a variety of 
clean metal surfaces and on MgB 2 (0001) were reproduced 

§6], |7|, g ||, H> H' l|L In a more sophisti- 
cated approach, Eguiluz et al. |64| obtained the correct 
asymptotic behaviour of the exchange-correlation poten- 
tial at large distances outside a jellium surface by solv- 
ing the so-called Sham-Schulter integral equation, which 
relates the exact exchange-correlation potential of DFT 
to the electron self-energy of many-body theory |55). A 
parametrization of this non-local potential was then used 
to obtain a Rydberg series of image states in the (100) 
and (111) surfaces of Al and Pd J66[ . 

Alternatively, the excitation energies and quasiparticle 
wavefunctions of unoccupied states can be calculated rig- 
orously from a one-particle Schrodinger equation contain- 
ing an energy-dependent non-local complex self-energy 
|p7[ . Simplified free-electron gas (FEG) models of the 
electron self-energy were employed by Echenique and col- 
laborators |38[ |6^, (70[ [7l] to calculate both binding en- 
ergies and a local effective potential. Full GW |7^] self- 
energy calculations of the binding energy and the local 
effective potential were reported in Refs. |73|, |74| for the 
lowest-lying image resonance at a jellium surface, show- 
ing that the effective potential nearly coincides with the 
exchange-correlation (xc) potential of DFT obtained in 
Ref. |Q. Furthermore, it was shown in Refs. (7^, Q 
that upon inclusion of long-range correlations into the 
xc potential |S3 the DFT eigenfunctions and eigenval- 
ues for a jellium surface are extremely good approxi- 
mations to their quasiparticle counterparts. More re- 
cently, full GW self-energy calculations for quasiparticle 
states at the (111) surface of real Al where reported |j75| , 
within 1.5 eV of the vacuum energy, which led to a first- 
principles evaluation of the image-plane position for this 
surface. The image plane for the many-electron system 
was found to be closer to the surface than that for the 
classical response to external fields of charges j7^| , and it 
was also found to be significantly modified by the atomic 
structure of the surface. 

The finite lifetime of excited states can be obtained 
from the knowledge of the imaginary part of the electron 



self-energy |7j . As the center of gravity of image states is 
located outside the solid, Echenique and Pendry [jll) as- 
sumed that the broadening of these states is dictated by 
their small overlap with the bulk region of the surface; 
for a constant value of the imaginary part of the bulk 
self-energy, they found that the image-state broadening 
is for all the levels in the Rydberg series smaller than 
the level spacing, thereby demonstrating the integrity of 
image states at metal surfaces. By combining TR-TPPE 
with the coherent excitation of several quantum states 
six image states were recently resolved on a (100) surface 
of Cu 1 37 , and their linewidths were found to approxi- 
mately increase with the quantum number as originaly 
predicted in Ref. |ll| . 

The first quantitative evaluation of image-state life- 
times was reported in Ref. j78|. In this calculation, 
the image-state wave functions were approximated by 
hydrogenic-like wave functions with no penetration into 
the solid and a simplified FEG model was used to ap- 
proximate the electron self-energy. In subsequent calcu- 
lations the penetration of the image-state wave function 
into the crystal was allowed |79, 80 1, and the role that 
the unoccupied part of the narrow Shockley surface state 
on the (111) surfaces of Cu and Ni plays in the decay 
of the n — 1 image state on these surfaces was inves- 
tigated by Gao and Lundqvist pi] . In this work, the 
image-state wave functions were also approximated by 
hydrogenic-like wave functions with no penetration into 
the solid, a simplified parametrized form was used for 
the Shockley surface-state wave function, and screening 
effects were neglected altogether. A GW calculation of 
the imaginary part of the electron-self energy near a jel- 
lium surface was also reported [ p2[ , showing the key role 
that a full evaluation of this quantity may play in the 
description of surface-state lifetimes. 

The first self-consistent many-body calculations of 
image-state lifetimes on noble and simple metals were 
reported only recently |8^, |84| , and good agreement with 
experimentally determined decay times ]3^, [38|, ^9) was 
found. In these calculations all wave functions and ener- 
gies were obtained by solving a one-particle Schrodinger 
equation with a realistic one-dimensional model poten- 
tial, the potential variation in the plane parallel to the 
surface was considered through the introduction of an 
effective mass, and the electron self-energy was evalu- 
ated in the GW approximation. Later, self-consistent 
calculations of the key role that the partially unoccupied 
Schokley surface state plays in the decay of image states 
on Cu(lll) were carried out j85|, and the inclusion of 
short-range xc effects was also investigated p6[ ]. It was 
demonstrated that although the presence of short-range 
exchange and correlation between screening electrons sig- 
nificantly enhances the decay probability of image states, 
this enhancement happens to be more than compensated 
by the large reduction in the decay rate produced by the 
presence of a xc hole around the image-state electron it- 
self. 

Recently, photohole lifetimes of Shockley surface states 
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were investigated in a variety of metal surfaces with the 
use of high-resolution angle-resolved photocmission spec- 
troscopy HI, HI H H © Hi- Scanning tunneling 
spectroscopy was also used to determine the lifetime of 
excited holes at the edge of the partially occupied Shock- 
ley surface-state band on the (111) face of noble metals 
p3[ , and to measure the lifetime of Shockley surface- 
state and surface-resonance electrons as a function of 
their energy Q. GW calculations |95|, |9(], |97j demon- 
strated that the decay of Shockley surface-state holes is 
dominated by two-dimensional electron-electron interac- 
tions screened by the underlying three-dimensional elec- 
tron system, and electron-electron interactions were then 
combined with the electron-phonon coupling to find good 
agreement with the experiment |ps| . These theoreti- 
cal investigations were extended to study the dynamics 
of surface-state and surface-resonance unnocupied elec- 
tronic states in Cu(lll) [ p8| , showing that, contrary to 
the case of surface-state holes, relatively major contri- 
butions to the e-e interaction of surface-state electrons 
above the Fermi level come from the underlying bulk elec- 
trons, and thereby giving an interpretation to the mea- 
surements reported in Ref. |94[ . 

Unless otherwise is stated, atomic units are used 
throughout this paper, i.e., e 2 = h = m e — 1. The atomic 
unit of length is the Bohr radius, ag = h 2 /m 2 = 0.529A, 
the atomic unit of energy is the Hartree, 1 Hartree = 
e 2 /ao = 27.2 eV, and the atomic unit of velocity is the 
Bohr velocity, vo — ac — 2.19 x 10 8 cms~ 1 , a and c be- 
ing the fine structure constant and the velocity of light, 
respectively. 



II. BASICS OF IMAGE STATES 

Image states are quantum states trapped in the image- 
potential well outside the surface of a material with nega- 
tive electron affinity. In the case of a hard-wall substrate 
that occupies the half-space z < and has a static di- 
electric constant e, the asymptotic form of the potential 
experienced by an electron in the half space z > is the 
classical image potential 



W 4z e + 1 



(1) 



are the eigenf unctions and eigenvalues of the one-electron 
Schrodinger equation 



ij){z) + V{z) =eiP(z). 



(4) 



Introduction of Eq. (||) into Eq. (^) gives rise to a Ryd- 
berg series of image-potential induced bound states (see 
Fig. 1) with 
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32 n 2 e+ 1 



n = l,2,. 



(0) 



where ip^ &en (z) represent the well-known wave func- 
tions of all possible s-like (I = 0) bound states of the 
hydrogen atom. Hence, the characteristic 'Bohr radius' 
of image states is expected to be 



4 e - 1 



a-B = — n — 
n e 
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do 



n=l,2,. 



(7) 



In the case of liquid helium, the static dielectric constant 
is very close to unity (e = 1.0572 for 4 He [Q), the bind- 
ing energy of image states is therefore of the order of a 
few tenths of meV (e = — 0.66/n 2 meV), much smaller 
than the actual barrier height Jl00[ , and the assumption 
of an infinitely repulsive substrate seems reasonable. 

In most metals there is no repulsive barrier for an elec- 
tron to return to the metal, but in a few metals there 
is a band gap near the vacuum level which may trap 
an electron in the quantum states of the image poten- 
tial (see Fig. 2). For metals the static dielectric con- 
stant is infinitely negative (e — > — oo), and image-state 
binding energies are of the order of a few tenths of eV 
(e = — 0.85/n 2 eV). As these energies are comparable 
with the typical width of the band gap, the condition 
of an infinitely repulsive substrate must be removed. A 
simple way to overcome this dificulty can be found in 
the multiple scattering approach developed in Ref. [pr| , 
which pictures surface states as standing waves of an elec- 
tron bouncing back and forth between the bulk crystal 
and the surface barrier. 



If the substrate is assumed to be infinitely repulsive, then 
the electron wave functions and energies are of the form 



1 



and 



*(r) = -=: V>(z)e ik " 
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E = e+-kl 



(2) 



(3) 



where k|| and rii are wave and position vectors parallel to 
the surface, A is the normalization area, and ^p(z) and e 



A. Multiple scattering approach 

Consider an electron wave ip+ propagating freely from 
the crystal termination at z = 0, which we choose to be 
located at one-half the atomic spacing beyond the last 
atomic layer, to the surface barrier at z = z- lm + z c {z- lm is 
the image-plane position), as shown in Fig. 3. This wave 
will be Bragg reflected both at the crystal [z = 0) and at 
the surface (z = z- lm + z c ) barrier. After both reflexions, 
the total reflectivity will be 
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r c n e 



i[0 B +2 K{z ira +z c )+4>a] 



(8) 
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where 4>b and 4>c denote the phase changes between inci- 
dent and reflected waves that are produced by the crystal 
and the surface barriers, respectively, and k is the per- 
pendicular component of the free-electron wave vector at 
< z < z- lm + z c . Summing the repeated scatterings gives 



1 - r B r c e 



i [</>c+2 K (zi m +z c )+4> b ] 



(9) 



Bound states correspond to the poles of the total ampli- 
tude dictated by Eq. (||), and occur when 



tb = rc = 1 



and 



4> B + 2 k (z im + z c ) + 4> c — 2nn , n = 0,1,2,. 



(10) 



(11) 



The phase change produced by any general potential 
barrier is obtained by matching the wave function ip(z) 
and its derivative ip{z) at the boundary of the barrier 
with the total (incident plus reflected) wave. One finds 



IK — L 



axg- 



i k + L ' 



(12) 



where L is the logarithmic derivative L = ip/ip of the 
wave function at the boundary of the barrier. 



1. Infinite crystal barrier 

In the limiting case of an infinite crystal barrier at 
z = and a pure image potential at z > 0, the crystal- 
induced phase is easily found to be <pc = 7r, k is infinite 
for finite energies, z- lm = z c = 0, and the cigenfunctions 
of the one-particle Schrodinger equation that are finite 
at infinite are the Whittaker functions 



4>(z) = W XA/2 (z/2X), 



where 



A = (-32 e) 



-1/2 



(13) 



(14) 



For A = n (n = 1,2,...), W A1 / 2 (0) = 0, L is infinite, 
and Eq. ( |l2| ) yields 4>b = (2n — 1) n, while for A ^= n 
(n = 1,2,...) Eq. fll2] ) yields 4>b = 2nir. Hence, in 
this limiting case ^ is a step-like function of e and the 
Bohr-like quantization condition of Eq. (|TT| ) is satisfied 
by the Rydberg series of bound states dictated by Eq. 
with e — > co, i. e. 



1 



32 n 2 



n = 1,2, 



(15) 



The n = condition of Eq. (O) is never satisfied 
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FIG. 3: Schematic representation of the full one-electron po- 
tential together with the image potential and image plane 
position Zim- 



2. Abrupt step barrier 

Alternatively, one may consider the limiting case of an 
abrupt step barrier of height Vq at z = z c . In this case, 



b B = 2 arctan j- [-e/(e - V )} 1/2 | 



(16) 



and the n — condition (0^ = —<Pc) can be satisfied, 
which yields the well-known crystal-induced surface state 
that exists in the presence of a Shockley-inverted gap. 
Image states are not present in this model. 



3. Quantum defect 

In the more general case of an image potential of the 
form 



V{z) 



1 



4(2 







(17) 



that joins the inner potential Vq at z = z c , both the 
crystal-induced n — Shocklcy state and the infinite 
Rydberg series of image states may exist. While the n = 
solution occurs at energies that are close to the Fermi 
level, the n / solutions occur at energies that lie near 



the vacuum level. For these energies, the variation of the 
crystal-induced phase 0c is small and the variation of the 
surface-barrier induced phase s can be descri bed b y the 



Wentzel-Kramers-Brillouin (WKB) expression[101 
4>b = (2/A 



1 7T. 



(18) 



Neglecting the energy dependence of 0c, substitution of 
Eq. @ into Eq. (y) yields 



32(n + a) 2 



1,2, 



where 



0c /tt) 



(19) 



(20) 



Here, the so-called quantum defect parameter a is re- 
lated to the bulk band structure through the crystal- 
induced phase 4>c- In the limiting case of an infinite 
crystal barrier, cf>c = ft and the quantum defect parame- 
ter is a = 0. 

In order to take explicit account of the variation of 
0c with energy and to describe the binding energies of 
both crystal-induced (n — 0) and image-potential in- 
duced (n 7^ 0) surface states, Smith combined the mul- 
tiple scattering description of these states with a NFE 
description of the bulk band structure and the use of 
only two bands 



III. TWO-BAND MODEL 

In the two-band approximation to the NFE band struc- 
ture of a solid, electron energies E(k) are the solutions 
of 



E(k) ~E°(k) 



-V* 



E(k)~E°(k-g) 



(21) 



where E°(k) = k 2 /2 is the free-electron energy and V s is 
the pseudopotential coefficient associated with the recip- 
rocal lattice vector g. At the k = g/2 point the degener- 
ate free-electron levels E°(k) and E°(—k) are separated 
by a band gap of magnitude 2|Vg|. Within this gap, solu- 
tions of Eq. (|l|) with real energy are still possible with 
a complex k z = p + i q corresponding to wave functions 
ip(z) that decay away from the surface into the bulk. 

In the case of a gap that is opened by potential Fourier 
components corresponding to g vectors that are normal 
to the surface, Eq. ( [H|) yields 



P = 9z/2, 



q 2 /2=J4(e-V Q )E°( g /2) + VZ 



(22) 

[e + E ( S /2)-V } , 
(23) 
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and 



= e q( - z 



" Zo) cos(p(z- z Q ) + 5). 



(24) 



Here, zq is the z coordinate of a surface atom, Vq is the 
crystal inner potential, and S represents a phase param- 
eter, 



sin 26 



(25) 



which in the presence of a Shockley-inverted band gap 
(V s > [102 1) varies from —ir/2 at the bottom of the 
gap to at the top of the gap. Matching at z — to a 
wave function of the form 



one finds, 



and 



ip(z) = e iKZ +r c e [ 



r c = 1 



K tan(0 c /2) = | tan 



(26) 



(27) 



(28) 



Hence, in the case of a Shockley-inverted band gap the 
phase <pc varies from at the bottom of the gap to 7r at 
the top, and the quantum defect of Eq. ( po[ ) varies from 
a = 1/2 at the bottom to a — at the top. 

In the more general case of a band gap that is opened 
by potential Fourier components corresponding to g vec- 
tors that are not necessarily normal to the surface, and 
at the points of high symmetry at the edge of the surface 
Brillouin zone where k x = g x /2 and k y = g y /2, there 
are two branches (j)^ and (j)^ depending on whether the 
wave function <j>(z) within the crystal is of even or odd 
symmetry with respect to surface atoms, 



k tan i 



6+/2) = |tan(|+* 



and 



k tan (cj> c /2) = -| cot (f + <*) • 



(29) 



(30) 



This shows that two surface states may occur in the same 
gap, and in cases where two such surface states are ob- 
served their energy separation is a measure of the surface 
corrugation potential, as pointed out by Smith [Q . 

By first introducing the experimentally determined 
band gap parameter V g and crystal inner potential Vb 
into either Eq. (||) or Eqs. (||) and Q, and then 
using the explicit energy dependence of <pB dictated by 
Eq. (Q, Smith and collaborators @ |9j solved Eq. 
( |ll| ) for the energies of the surface states of a variety 
of face-centered cubic metals. Both crystal-induced and 
image-potential induced surface state binding energies 
were found to be in reasonable agreement with photoe- 
mission and inverse photoemission measurements, and 
the systematics of surface-state occurrence in bulk band 



gaps on different crystals were examined in an elegant 
way. 

Instead of using the WKB expression for <pB [Eq. 
which is only appropriate for energies that lie very near 
the vacuum level, Ortuho and Echenique [ |50| considered 
the more accurate expression 



t>B = 2 tan 



1 J (y) cos(Att) + N (y) sin(ATr) 



Ky Ji (y) cos( A %) + Ni (y) sin(A n) 



(31) 



as obtained b y int roducing into Eq. (|12| ) the Wannier 
approximation ]103[| to the Whittaker functions, i.e., 



W x , 1/2 (y) « y [My) cos(A^) + iVi(y) sin(ATr)] . (32) 



Here, y = ^2 (z~ + z c ), and J (y), Ji(y), N Q (y) and 
N\(y) are Bessel and Neumann functions, respectively. 

By using Eq. (pi]), Ortuno and Echenique ||5(| solved 
Eq. ( Jll| ) for the energies of the surface states associated 
with the bulk L gap projected onto the (111) surfaces 
and the bulk X gap projected onto the (100) surfaces 
of face-centered cubic Cu, Ag, and Ni. It was demon- 
strated that this model correctly predicts the trends and 
systematics of both crystal induced and image-potential 
induced surface states, and was found that the experi- 
mental binding energies are reproduced for image-plane 
positions that are closer to the crystal surface than those 
predicted by Lang and Kohn|[76[, 



IV. ONE-DIMENSIONAL MODEL POTENTIAL 

At large distances from the surface, the one- 
dimensional potential V(z) experienced by an electron 
has the asymptotic form exhibited in Fig. 3, with the 
image plane located at z — z lm . However, on approach- 
ing the surface V(z) must merge continuously into the 
bulk crystal potential. Jones, Jennings, and Je psen pro- 
posed a saturated model potential of the form [104] 



V(z) 



1 +4c°( 2 - z '») 
1 



4 (Z - Zi m ) 



1 - e 



-A(z-z im ) 



Z < Zj, 



Z\ » T i ^ z , 



(33) 



where Vq is the inner potential, A -1 is a characteristic 
distance for the changeover between the inner potential 
and the image asymptotic form, and the parameters A 
and a are fixed by the requirement of continuity. By nu- 
merically integrating the Schrodinger equation with the 
model potential of Eq. (|33"1), Smith and collaborators |5l[] 
obtained the phase change 4>b and applied the two-band 
model described above to estimate the image-plane posi- 
tion of a variety of metal surfaces. They found that while 
the face dependence of the image-plane position is a di- 
rect result of the discrete lattice nature of a real material, 
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the response of the electron density is the fundamental 
quantity determining the image-plane position itself. 

Recently, a model potential which involves the use of 
the n — and n — 1 surface-state binding energies as well 
as the width and positio n of the energy gap as fitted pa- 
rameters was proposed [105]. This one-dimensional film 
model potential has the form 



' Aio + A\ cos(27rz/a s ) z < D 

A 2Q + A 2 cos [ft{z-D)\ D<z< zi 

A 3 e- a ^- z ^, *i<z<z im 
1 



V(z) = I 



I 4(z-z im ) 



1-e 



■\(z—z lm ) 



%\ m ^ Z . 



(34) 



where the origin of coordinates has been taken to be lo- 
cated at the middle plane of the film, D is the half-width 
of the film, and a s is the interlayer spacing. This model 
potential has ten parameters, A 2 $, A3, a, z\, A, z- lm , 
Aio, Ai, A2, and ft, but only four of them are indepen- 
dent. A20, A3, a, z\, A, and z- lm are determined from the 
requirement of continuity of the potential and its first 
derivative everywhere in space, and the remaining four 
parameters are chosen as adjustable parameters. The pa- 
rameters A\ and A\q reproduce the width and position 
of the energy gap, while A 2 and ft are chosen to repro- 
duce the experimental or first-principles n — and n = 1 
surface-state binding energies at the f (k|| = 0) point. 
Although band-structure calculations that are based on 
the model potential of Eq. ( |34| ) yield work functions 
that are within 0.1 eV of the measured values, for a more 
precise description of the unoccupied states below the 
vacuum level the use of experimental work functions is 
recommended. 

The model potential of Eq. ([34]) was constructed in 
Ref. 1 59 for 14 simple and noble metal surfaces. By using 
this model potential both wave functions and binding 
energies of image states were investigated, and the issue 
of the image-plane position was also addressed. Image 
planes were found to be typically closer to the crystal 
surface than those predicted by Lang and Kohn U& for 
a jellium surface, and they were found to be, in the case 
of noble metals, in agreement with those reported in Ref. 
filf by Smith and collaborators. Good agreement was 
also observed with the image-plane positions reported in 
Ref. 1 66 from first-principles band-structure calculations 
of the (111) and (100) surfaces of Pd. 

Key quantities in the description of the lifetime of sur- 
face states are surface-state wave functions and charge 
densities. As shown in Ref. |5^, the model poten- 
tial of Eq. (j34| ) yields charge densities of both n = 
and n = 1 surface states that are in excellent agreement 
with those obtained with the use of first-principles calcu- 
lations. 



V. LIFETIMES 

Accurate im age- s tate lifetimes were reported in Refs. 
|[S3"| , |84| , |85| , [86 , 106, 107 1, as obtained from self-consistent 
many-body calculations involving the use of the one- 
dimensional model potential of Eq. (|34j). 

In the framework of Green-function theory, one iden- 
tifies the inverse quasiparticle lifetime as follows 

r7 1 =-2ImS s , (35) 

where E s is th e qu asiparticle energy. On the energy-shell 
one finds I IPS] 



r- 1 = 



2 J dr J dr'*:(r)lml](r,r';£ s )* s (r), (36) 



where ^s(r) and E s are taken to be the eigenfunctions 
and eigenvalues of a suitable Hermitian single-particle 
Hamiltonian and £(r, r';E s ) is the so-called self-energy of 
the quasiparticle, which is a nonlocal, energy-dependent, 
non- Hermitian operator accounting for all xc effects be- 
yond the Hartree approximation. Assuming translational 
invariance in the plane of the surface, i.e., taking ^(r) 
and E s to be of the form dictated by Eqs. @ and ||), 
Eq. <M) yields 



dz I dz' t/}*(z)ImE(z,z';'k\\,E s )4> s (z), 



(37) 



where ijj s (z) and e s represent quasiparticle wave func- 
tions and energies describing motion normal to the sur- 
face. In usual practice, these wave functions and energies 
are taken to be either the eigenfunctions and eigenvalues 
of a one-dimensional LDA Kohn-Sham Hamiltonian of 
DFT or the solutions of a one-dimensional Schrodinger 
equation with the model potential of Eq. (F34|) . 



A. GW approximation 

The exact self-energy E(z, z'\ k|| , e s ) can be obtained, 
in principle, from an iterative solution of Hedin's equa- 
tions p7| in combination with the Dyson equation [ 109 1 . 
However, to obtain explicit results one usually resorts to 
an expansion in powers of the time-ordered screened in- 
teraction W(z, z'; kii, £■). The leading order term of this 
expansion is the so-called GW approximation: 



f°° de 

£(z,z';k||, Es )= / — e-^G(z,2';k||, £s - £ ) 
J-00 2tt 

x W(z,z!]\L\\e), (38) 

which can also be obtained as the first iteration of 
Hedin's equations by simply neglecting vertex correc- 
tions. G(z, z'; k|| , e) represents the one-particle Green 
function, rj is a positive infinitesimal, and the screened 
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interaction can be expressed in terms of the density- 
response function xi z i z ')^\\i e )i as follows 

W(z, z';ky,£) = v(z — 2',k||) + / dz x J dz 2 
xv(z - zi,k||) x(zi,z 2 ;'k\\,e) v(z 2 - z). (39) 

Most current GW calculations simply replace the exact 
one-particle Green function entering Eq. ( |3^ ) by the non- 
interacting Green function 



G°(z,z';k||,e)=^ 



e - Ef + i?7Sgn(e - E F ) 



(40) 



where ipf(z) and £/ are taken to be a complete set of ei- 
ther eigenfunctions and eigenvalues of a one-dimensional 
LDA Kohn-Sham Hamiltonian or solutions of a one- 
dimensional Schrodinger equation with the model poten- 
tial of Eq. (||). By introducing Eq. @ into Eq. (§£ 
and then Eq. (|38j) into Eq. (p7|), one finds 



= - 2 £ 



dz' 



xlmW(z, z'; q„, u)<t> f (41) 



On the same level of approximation and neglecting all 
vertex corrections, the screened interaction entering Eq. 
(fll]) is usually obtained from Eq. (|3^) with the density- 
response function evaluated in the random-phase approx- 
imation (RPA): 

X0> z';k||,e) = x°(z,z';k||,e) + / dzi / dz 2 

x x°(z, z '-, k || i e) v{zi - z 2 ; ky , e) x(z 2 , z'; k|, , e),(42) 

X°(z, z'; k|| , e) being the noninteracting density-response 
function 



X °(z,z';k h e) 



-2i 



de' G°(z,z';k||,e') 
xG°(z,z';k| h e + £')- (43) 



B. GWT approximation 

Short-range xc effects, which are absent in Eqs. (El]) 
and (|4^), can be included i n the fram ework of the so- 
called GWT approximation 110 , 111]. In this approxi- 
mation, the electron self-energy is of the GW form, i.e., 
it is given by Eq. (|38|), but with an effective screened 
interaction 



W(z,z';k\\,s) = v(z - z';k||) + / dzi / dz 2 

x [v(z - zi;k||) + / xc (zi,z 2 ;k||,e)] 
xx(2i,z 2 ;k||,£)u(z 2 - z';k||), (44) 



the time-ordered density-response function 
X(zi,Z2','k\\,e) now being 



X(z,z';k||,e) = x°(^,^';k||,£) + / dz x / dz 2 

xx°(z,z';k||,£) [w(z - z';k||) + / ac (zi, 2 2 ; ky, £)] 
x x (z 2 ,z';k h e). (45) 

The kernel f xc {z\, z 2', k||, e) entering Eqs. (^3) and fl45| ) 
accounts for the reduction in the electron-electron in- 
teraction due to the existence of short-range xc effects 
associated to the excited quasiparticle and to screening 
electrons, respectively. In the so-called time -dependent 
local density approximation (TDLDA) |112| or, equiv- 
alently, adiabatic local-density approximation (ALDA), 
the exact xc kernel is replaced by 



fxc ( z > z > k|h e ) 



dn 2 



5(z - z'), 



i(z) 



(46) 



where £"™^ (n) is the xc energy per particle of a uniform 
electron gas of density n, and n(z) is the actual electron 
density at point z. 



C. Effective mass 

The dispersion E(kn ) of both bulk and surface states 
has been determined experimentally with the use of in- 
verse photoemission techniques at off-normal emission 
filf , showing that it is of the form dictated by Eq. (||) 
with k||/2 replaced by kii/(2m), m representing an ef- 
fective mass. Surface corrugation effects on the effective 
mass of image states have been found to be negligible, 
and many-body effects are found to enhance the effective 
mass by no more than « 2%. Nevertheless, the effective 
mass of the n = crystal induced surface state and that 
of unoccupied bulk states considerably deviate from the 
free-electron mass, and the impact of this deviation on 
the lifetime of both crystal-induced and image-potential 
induced surface states is found to be important. On 
the one hand, there is the effect of the decrease of the 
available phase space, which is easily found to scale as 
■JmJ, nif being the effective mass of the various avail- 
able final states. On the other hand, as the effective 
mass decreases, the decay from the image state occurs, 
for a given energy transfer, through smaller parallel mo- 
mentum transfers, which may result in either enlarged or 
diminshed screened interactions, depending on the mag- 
nitude of momentum and energy transfers. 



D. Results 

In Table 1 we compare theoretical results with 2PPE 
and TR-2PPE data for image states at the T point on 
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TABLE I: 
meV. 



Linewidth (inverse lifetime) of image states, in 









2PPE 


TR2PPE 


Theory 


Cu(100) 


ii= 


= 1 


28 ± 6° 


16.5 ±3/2,3'' 


22 e ; 17 1 




ii= 


=2 




5.5 ± 0.8/0.6,3 6 


5 e 




ii= 


=3 




2.20 ± 0.16/0.14,3 6 


1.8" 


Cu(lll) 




;1 


85 ± 10 a 


38 + 14/9 9 - 30 h 


38 c ;29 / 


Agl^lUUJ 


ii= 


- 1 

- 1 


91 4- zl a 


1 O _|_ 1 o*> 

±z zrz ±,o 


on/c 

00 




ii= 


-0 


o. i m u.^t 


41 _i_ n /n 9 o& 
^.1 zrz u.o/u.Zr,o 






ii= 


_q 
-o 




l.OO ztz U.Uo 




A „/l 1 1 \ 

Ag(lll) 


ii= 


-1 


40 ± 1U ;oo 


00 1 t n in tn 
AZ dz 1U/D 




a nl ^ i i i ^ 
Au^lll ) 


ii= 


i 

-1 


LOU zh 4U 






lNl(lUO) 


ii= 


i 

= 1 










n= 


= 1 


84 ± 1(J 






Co(0001) 


n= 


= 1 


95 ± 10 a 






Fe(llO) 


n= 


= 1 


130 ± 30 a 






Pt(lll) 


n= 


= 1 




25 ± 10/5" 


29™ 




n= 


=2 




11 ± 1/2™ 


9" 


Pd(lll) 


n= 


= 1 


70 ± 8 a ;32° 


27° 


30° 




n= 


=2 






7° 


Li(llO) 


n= 


= 1 






37 p 




n= 


=2 






15 p 



"From 
'From 
c From 
c From 
^From 
'""From 
m From 
"From 
"From 
p From 



Ref. 31 

Ref. 22. 

Ref. S 

Ref. S 



Ref. 
Ref. 
Ref. 
Ref. 
Ref. 
Ref. 



so 



various metal surfaces. While the calculations reported 
in Ref. J8j| were carried out within the GW approxima- 
tion with use of the free-electron mass, realistic effective 
masses were i ntro d uced in the calculations reported in 
|l(J7. 1 13 1 . GWT calculations were only re- 



Mi 



106 



Refs. 

ported in Ref. |86| for the n = 1 image-state lifetime on the 
(100) and (111) surfaces of Cu. It was found that differ- 
ences between the calculations reported in Refs. js^. |8(J 
for these surfaces were mainly due to the impact of the 
actual effective mass on the decay mechanism, and was 
shown that GW calculations produce decay rates that are 
very close to GWT calculations. Both GW and GWF cal- 
culations produce image-state lifetimes that are in good 
agreement with TR-2PPE data for Cu(100), Cu(lll), 
Pt(lll), and Pd(lll). Even for Ag surfaces, where the 
polarization of rf-electrons is important for the descrip- 
tion of the surface and bulk plasmons and which has not 
been included in the calculation, GW calculations lead 
to a reasonable value for the image-state lifetime. Calcu- 
lations that include the effect of d-electron polarization 
are now in progress [114 1. 
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